ITPSB 95-14 
HEP-TH 9509004 



Analyticity and Integrability 
in the 

S : Chiral Potts Model 



00 



X 
J3 



Barry M. McCoy and William P. Orrick I 



Stony Brook, NY 11794-3840 



^ ! Institute for Theoretical Physics 

^ ■ State University of New York 

o ■ 

O 

in 



9/95 



^ mccoy@max.physics.sunysb.edu 
^ worrick@insti.physics.sunysb.edu 



Abstract 

We study the perturbation theory for the general non-integrable chiral Potts model 
depending on two chiral angles and a strength parameter and show how the analyticity of 
the ground state energy and correlation functions dramatically increases when the angles 
and the strength parameter satisfy the integrability condition. We further specialize to 
the superintegrable case and verify that a sum rule is obeyed. 

1. Introduction 

The connection between integrability and analyticity has been recognized and ex- 
ploited for many years in the theory of solvable models in statistical mechanics and in par- 
ticular in the method of solution based on the inversion relation where assumptions 
on analyticity are a vital ingredient of the solution. Conversely, whenever an integrable 
model can be solved by other means the analyticity demanded by the inversion relation is 
always found to exist. 

On the other hand these analyticity statements have never been made particularly 
precise and, indeed there are recent computations in certain four dimensional field the- 
ories 0-0 where results may be obtained because functions also have dramatically ex- 
panded analyticity properties over what is generically expected but no analogue of the two 
dimensional integrability conditions is known. Consequently it is of interest to sharpen 
our understanding of the connection between integrability and analyticity. 

In this note we examine these questions for the state chiral Potts spin chain defined 

by 

The matrices Xj , Zj are defined by 

Xj = In®---^^---^In , Zj = In ®---^Z^---^In , (1-2) 

site J site j 

where In is the N x N identity matrix, the elements of the N x N matrices X and Z are 

Xl,m = (mod A^) , Zi^rn = (1-3) 

and uj — e^'^'^l^ , When the angles and the strength parameter k satisfy the condition 

cos (f) = k cos (j) (1-4) 
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this spin chain Hamiltonian is derivable from the state chiral Potts model whose transfer 
matrix Tp g satisfies the integrability condition of commuting transfer matrices 

[Tp,,.Tp,,']=Q (1.5) 
where p and q specify points on the Riemann surface given by the intersection of two 
Per mat cylinders 

+ A;6^ = A;'c^^, A;a^ + 6^ = A;'c^ (1.6) 

with /c^ + /c'^ = 1. The Hamiltonian ( p. . 1|) is obtained from Tp ^ in the limit g — > p as [|] 

Tp,g = / + w(ifcp + const) +0(^2) (1.7) 
where u measures the deviation of q from p and 

bpdp bpCp 

The free energy of the statistical mechanical system and the ground state energy of 
the spin chain have been computed in the general integrable case [§] and in the special 
superintegrable case of^ = ^ = 7r/2 where there is a dramatic simplification H ||10|| - 
||11|| . Our study here is to relate the general non-integrable case to the integrable and 
superintegrable cases. 

Por the non-integrable case the only analytic tool available is perturbation theory in 
the variable k. This expansion can be done both for /c ~ and A; ~ oo. In this paper we 
study only the first case and consider the expansion for the state spin chain of the ground 
state energy and correlation functions. We observe that as the order increases, new singu- 
larities in the e*"^ plane continually arise and we expect that the number of singularities will 
become infinite as the order goes to infinity. The integrability condition ( |1.4| ), however, 
involves k and causes extensive cancelations among orders. We find that all but a finite 
number of singularities, which are already present in the lowest orders, vanish resulting in 
a dramatic increase in analyticity. This increase in analyticity is intimately related to the 
integrability of the model. 

In Section 2 we present the perturbation expansion of the ground state energy for 
general chiral angles cj) and (j) for = 3-6 and then exhibit the increase in analyticity that 
occurs when the integrability condition ( p..4|) is imposed. In Section 3 the same procedure 
is applied to the correlation functions for = 3. In addition, the perturbation expansion 
for the superintegrable correlation functions is given for arbitrary A^. In Section 4 a sum 
rule for the nearest neighbor superintegrable correlation functions is derived and used to 
verify the expansion in Section 3. Also, a Puchsian equation for the superintegrable ground 
state energy is derived. In Section 5 we conclude with some questions for further study 
and some general remarks on the connection between integrability and analyticity. 
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2. Perturbation Expansion of the Ground State Energy 

In this section we consider the expansion for the state chiral Potts spin chain with 
general chiral angles, (f) and (f), of the ground state energy 



e(^) = hm eSVl = $:fc-e(^). (2.1) 



L — >oo 

n=0 



We have implemented Rayleigh-Schroedinger perturbation theory [|T^ on the computer 



using a program written in C, Mathematica, Maple and Form and studied the following 
cases: 

A^=3 to order 9 N=A to order 8 A^=5 to order 5 A^=6 to order 6. 

Our results for these 4 cases are given in Tables 1-4 where we use the following notation 

C3 = cos(0/3), C4 = V2cos((^/2) 

C5 = ^^4^ = v^Vl + 5-V2 cos((/>/5), Ce = v^cos((/>/3). 
sm(7r/5) 



In the case of = 3 there have been previous perturbation studies both for small [|T3 



14] and large [14|-|[T1] k and Table 1 agrees with the small k results of |14] which were 



obtained to order 7. Many of the details of the perturbation expansion are discussed in 
these papers. 

These tables demonstrate an important feature of the generic perturbation theory; 
namely that as the order of perturbation theory increases so do the number of places in 
the e*"^ plane where the perturbation expansion develops singularities. Thus, for example, 
we see from Table 1 that 62^'* and e"^^ are singular only when C3 = but that 64"^^ Cg'^'', 
ef \ and e7'^''and have additional singularities at C3 = ±1/2 and ef^ and e^^ have further 
singularities at 28C3 — 3 = 0. A similar increase in number of singularities is seen in 
Table 2. Indeed this property of the number of singularities increasing with the order 
has been explicitly seen in a variety of previous perturbation computations including the 
Ising model in a magnetic field [|l| and the zero field Ising magnetic susceptibility fl^. In 
the latter case is it well understood that the number of singularities is connected to the 
number of quasiparticles that contribute to the given order of perturbation theory [jlS 



||19|| and that as the order of perturbation goes to infinity an infinite number of singularities 
will appear. It is this ever increasing number of singularities in the perturbation theory 
which in the case of the Ising model in a magnetic field restricts the analyticity to such an 
extent [|T1 that a closed form solution has never been found. 



However there is a dramatic change in this pattern of singularities when the parameters 



(/), ^, and k are related by the integrability condition (|1.4|) . This is possible because this 
integrability condition (|1.4|) involves the perturbation parameter k and thus allows different 
orders of perturbation theory to be combined. 

To appreciate the cancelations between the orders of perturbation theory that result 
we consider = 3 in detail. It is easy to see when ( |1.4| ) holds that if we use 



cos ( 



cos(0/3) 



i</)/3 



(2.3) 



then e^"^ and e'^' are combined as 



,(3) 



A.2p(3) , l3 (3) 



_^2 3 + 4C| 



973^3 



(2.4) 



which is the /c^ term on the integrable manifold as given in Table 5. 

(3) (3) 

A more important cancelation, however, takes place between the terms 64 , and 
e^^ to give the k"^ term in Table 5. Here the terms in cos (p in the k^ terms and the terms 
in cos^ (f) in the k^ term which have the denominators ±1 + 2C3 combine with the terms 
in the /c^ term with the same denominators to cancel the denominators out completely. 
In addition the terms with C^'^ in 65'^'' and with C^^ in e^^^ are reduced to C^^ by use 
of ( |2.3| ) and thus the complete k"^ term in Table 5 is obtained from three terms in Table 1. 

(3) ^ (3) 

Similarly, the terms in Cg which have no factor of cos(/) combine with the terms in Cy 
with cos the terms in eg with cos^ (p and the terms in e^^ with cos"^ (j) to cancel all 
denominators containing powers of ±1 + 2C3. Thus we obtain the term of order k^ in 
Table 3. 

For = 4 we find from Table 2 that the terms of orders 2, 3 and 4 are needed to get 
the /c^ term on the integrable manifold and that to get the /c^ term in the integrable case 
the terms in the generic case to order 8 are needed. In general to get the order k'^'^ term 
in the integrable case the terms to order Nq of the generic case are needed. To show this 
we first note that e*^^-* is an even function of (p. Secondly, the order r term in the generic 
case is a sum of terms of the form 



Cn,i{(t>)k'' exp 




1=1 







(2.5) 
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where |0)o is the ground state when k = 0. The product clearly vanishes if Yl^i=i '^i 
an integer multiple of A^. Thus the k'^ term can be written 

c» oo 

e(^) = ^k''Y^ Cp{(f>) cos{p4>) = XI ^'p^'f'^ ^^^'^ ^ (2.6) 

r=0 p r=0 p 

where the sum over p is taken over p even if r is even and p odd if r is odd. When the 
integrability condition ( [T^) is applied we obtain 

g(A^) =J2Y1 k^'-^Cpicj)) cosP (t> (2.7) 

r=0 p 

which contains only even powers of k. Since the maximum allowed value for the ni in (|2.5|) is 
— 1, the maximum allowed value of p for a given r is Pmax{i^) = 2 '^^"^ — r. The 
only terms in the generic case contributing to the /c^"^ term in the integrable case are those 
terms for which r — Pnia.x{r) < 2q which is equivalent to the condition — [— < Q or 
r < qN as desired. 

The results of the specialization to the integrable manifold are given in Table 5 for 
the cases = 3-6. Thus we see in the cases considered that all singularities have can- 
celed out except certain of the original singularities in 62"^^ . Consequently the analyticity 
of the ground state energy has increased to the maximum extent possible. This remark- 
able extension in analyticity is a consequence of the defining commutation relation of 
integrability ( p.. 51 ) and the relation ( |1.8| ) which the chiral angles have with the Riemann 
surface (|1.6|). 

These integrable results are now easily specialized to the superintegrable case (p = n/2 
and we obtain 

e(3) = -2 - - Afc4 _ 20^6 _ o(fc8) 



32 35 37 
5 2 25 
23 29 



,(5) = _4 _ 1^2 _ 0(^4) 

5 



e^'^ = -5 - - 0{k') 
36 



which agrees with the previous computation of |]T^. We note that in principle we should 
compare the integrable results of Table 5 with the exact result of but no systematic 
expansion of the exact result seems to be in the literature. 
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3. Perturbation Expansion of the Correlation Functions 

The cancelation of singularities and resulting expansion of analyticity seen in the 
ground state is, of, course, expected to occur in all the correlation functions as well. 
However, the presentation of the results of the perturbation theory in the generic case 
is more cumbersome and consequently we will restrict ourselves here to illustrating the 
phenomena for the case = 3. 

The perturbation expansion of the correlation function (ZoZj.) for generic values of (j) 
and d) was studied in [Ol-IITll to fourth order. In order to exhibit the desired cancelations 



of the poles at C3 = ±1/2, which occurred in the ground state energy by combining 
the fourth, fifth and sixth orders, we must extend the computation of the correlation to 
sixth order also. This we have done and the results are given in Table 6. The method 
of perturbation is discussed in some detail in [0|-[n]. The only additional remark we 



wish to make is that in contrast to more standard perturbation expansions we have been 
unable to find an efficient way to perform the subtraction of the disconnected diagrams 
to obtain a general formalism that deals only with order one objects in the L 00 limit 
and are forced to perform an explicit subtraction of L dependent terms. This seems to 
be a consequence of doing perturbation theory on a lattice with a variable that satisfies a 
cyclicity condition Xj^ = 1 and, because of the connection of the cyclicity condition with 
the fractional statistics of the model, it deserves further study. 

To study the specialization of the generic result of Table 6 to the integrable 
case ( p.. 41 ) we consider first the long range order term = liuir^ooiZoZ}.) which in 
the perturbation expansion is the coefficient of 1 — 6n,o- From the second and third order 
terms it is easy to see that the cancelation of poles is more complete than in the ground 
state energy and that to second order on the integrable manifold we have the simplification 
to — 2/c^/9 where the dependence on the chiral angle (p has completely vanished and thus 
the order parameter on the integrable manifold is identical with the order parameter on 
the superintegrable special case. This complete cancelation occurs in all orders and is a 
well known consequence [^-|]21[] of the integrability commutation relation ( |1.5| ). These 



perturbation computations for the long range order have been carried to high order in k 
for iV = 3 in [^T|-[^^ and have been extended to arbitrary in []T^. These expansions 



have lead to the remarkably simple conjecture for the general case that 

M2 = hm {{ZoZly) = (1 - A;2)"(^--)/^'. (3.1) 

~ +00 



To see the first dependence on the variable (f) in the correlation function on the inte- 
grable manifold we must go to fourth order which is obtained from the fourth, fifth and 
sixth orders of of the generic result of Table 6. (Note that to obtain the /c^'^ term in the 
integrable case we need terms up to k^'^ in the generic case just as we did for the ground 
state energy.) In fourth order the only deviation of {ZoZj.) from occurs for r = 0, 1 and 
since the case r = is trivial we consider only r = 1 and find that the poles at C3 = ±1/2 
cancel. Thus the result is obtained on the integrable manifold ( p..4| ) of 

7 1 



{ZoZl)-M! = k^{- 



(3.2) 



162 108cos2(^/3 ■ 54 
This result is to be compared with the fourth order result on the superintegrable 
manifold which may be obtained for arbitrary A^. From a straightforward perturbation 
expansion we may reduce the correlation to a set of trigonometric sums 



{(Zozir) 
1 



1 



+ 



16A^4 
1 



N- 



n 
N 

N-l N-1 

EE 

j = l k=l 
-1 



1 



n 
N 

2(1 - u""^ 



) + 2(l 



-nk 



)-(l 



UJ 



n{j-k) 



) 



8N^ 



E 



p=i 



1 



cos 



2nnp\ 



sin nj/N sin nk/N 

N-l N-l 



N 



EE 

/ ^ — 1 1 — 1 



j=i fe=i 



71 j nip 
sin — sin 



Txk T^ip — k) 



1 -1 



1 



N 



I (1 — cos2nnp/N) 



(3.3) 



4Ar4 



sin Txp/N 



+ 



1 



4iV4 



N-l 

E 

p=i 

N-l 



1 — COS 



2nnp 



N 



N-l N-l r 



sm 



N 



16A^4 



EE 

J = l k=l 



UJ 



EE 

nU-k)^ + (1 - a;-"J) + (1 - Uj'''' 



nk n(j — k) np 

sin — — sin — sin — — 

N N N 



-\ -1 



sin nj/N sin nk/N 



k^ + 0{k^ 



The trigonometric sums have been evaluated in ref. and thus the final result is 



{{ZoZir)-M^ 



— I 



n 
N 



1 - 



N-l 



n 
N 



+ 



1 



n 



4iV2 AT 



n 
N 



4iV4 



J2f{N,n,p)ts.n^jk^ + 0{k' 



(3.4) 



p=i 



where 

/(Ar,n,p) 
Thus if we set = 3 in 



n'^{N-2p) ifl <n<p< N/2 

f{N,p,n) in general 

—f{N, N — n^p) in general. 

I) we obtain (g^ with </> = f . 



(3.5) 



7 



4. Sum Rule for Superintegrable Correlation Functions 

The superintegrable nearest neighbor correlation functions {{ZqZI)'^) obey a sum rule 
which provides a check of the perturbation expansion ( |3.4D . Here we derive several forms 
for the sum rule. 



We start by applying the Feynman-Hellmann formula p4 

dH 



c>e(^) 1 , 



dk 



to (|1.1|) to obtain 



where 



N-l 



5](0|a^(ZoZ;r|0) = -e(^)+fc 



n=l 



dk 



(4.1) 



(4.2) 



,i{2n-N)<P/N 



ar 



sm{'Kn / N) 

and the ground state energy per site of the A^-state superintegrable chiral Potts chain is 



(4.3) 



N-l 



AN), 



1 / 



4k 



With the aid of |2|] (p. 102 (21)) 



(4.4) 



a 



F'{a, b;c;z) = - [F(a + 1, b; c; z) — -F(a, 6; c; z) 



(4.5) 



(O) and (O) yield a sum rule 



N-l 



N-l 



j2{oMZozir\o) = j2 



n=l 



n=l 



4k 



2' N' ' (1 + A;)2 



+ 2(1 + ^)^ 



1 n 



4k 



2' N' ' (1 + A;)2 



We may reexpress the sum rule in terms of the functions F{ 
F(— -S^, -S^; 1, A;2) by using the two identities 



TV' Af ' 



(4.6) 



1, k'^) and 



0; 1; 



l + kf^F{b, b;l;k^) 



{l + k)\ 

(-^, 1; (T^) =(1 + ^)''"' (4(1 - ^Mb, -b; l; fc^) 

+ [(46-3) -fc](l-A;)F(6,6;l;A;2)}. 



(4.7) 



The first of tfiese is just eq. (36), p. 113 of with a = b. The second can be obtained 
from 



=(26-1) 1 



4k 



{1 + ky 



, 1 , 4k 



1 4k 
+ 2(l-5)F(-,6-l;l;^^ 



(eq. (37), p. 103 of ^ with a = c = 1 and z = 4A;/(1 + A;)^) by using the first identity 
and eq. (33), p. 103 of [^S] with a = 6— 1, c=l and z — k'^. The sum rule now becomes 

N-l N-1 



iV— 1 iV— lr/„ s 

(OMZoZlriO) = ( f - 1 (1 - - k)^-'F (^, ^; 1; k 

n=l n=l L V / 



+ 2 1- 



n 
N 



(4.9) 



The two sides of the sum rule are even functions of k. This is made manifest by 



adding (|4.9|) with n N — n to itself and using eq. (2), p. 105 and eq. (31), p. 103 



of to obtain 



N-l 



El ■TT(2n — N) x ■7r(2n — JV) + 

-. ^(0|e^^^(ZoZ )- + e-^^^(ZoZ )^-|0) 
sm Tin IS 

n=l ' 



N-l 



5: {2(1-^) + + 



(4.10) 



n=l 



+ (|-i)(i-fc2) {i + k)^-' + {i-ky 



Expanding the right hand side of ( [4.10| ) to fourth order and using ( pi. 4] ) in the left 
hand side the k'^ term in the sum rule gives 



N-l N-l 

J2 fiN,j,k)cot ^^cot 
j=i k=i 



nj Tik _ N{N^ ~1){N'^ -4) 
iV ~ 90 



(4.11) 



This identity is a special case of the more general identity, which is independently proven 
in Appendix A 

N-l N-l . , 

s^^) = E E ^) ^ 77 

(4.12) 



j=i k=i 



2{N -1){N -2) 



g(iv,o)-4 Y {N~2j)Q[n,^ 

i=i 
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where 

(P{N,j,k) iil<j<k<N/2 
/Q(iV,i, k) = I /Q(iV, kj) in general (4.13) 

[ -/q(A^, N - j, k) in general 

and 

PiN.J. k) = Q (at, ^) + Q (at, ^) - 2Q (iV, A) (4.14) 
for any function Q{N,x) satisfying 

Q{N,l-x) = -Q{N,x). (4.15) 

We note that ( [4.12| ) is a member of a broader class of identities, described in Appendix B. 

We conclude this section by deriving an order 2N — 2 Fuchsian differential equation for 
the superintegrable ground state energy ( [4.4| ) with three singularities located at 0, 1 and 
oo. Equation (O) is expressed in terms of hypergeometric functions which are solutions 



to the second order Fuchsian equation with singularities as above 

I^(ate)i^(a,6;c;z) = (4.16) 

where 

(P d 

D(a,b-c-z) = z{l-z)-^ + [c-{a + b+ l)z] — - ah. (4.17) 

One may readily verify that this differential operator satisfies 

D{a+lM;c;z)D(^a,h2;c;z) = ^(a+1,62 ;c;2)^(a,f>i ;c;z) ■ (4-18) 

We also define the differential operator 

DU.;,) = (1 + ^)^(,,,;_^) (1 + k)-' (4-19) 

which has solution (1 + k)F ^a, 6; c; -^^^ j and which satisfies an equation precisely anal- 
ogous to ( p^) . Now let 

Dn=D'( 2N-5 N-1 .-..,\D'( 2N~7 N~2 ■ ■ ■ D'( l (4.20) 

By repeated application of ( [1.18| ) we see that 
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where (61, ... , bN-i) is any permutation of (1, . . . , — 1). In particular, bi may take any 
value from 1 to — 1 from which it immediately follows that 

DNe^^\k)=0. (4.22) 



5. Discussion 

No finite order of perturbation theory can ever stand as the final solution of any 
problem. However, such perturbation studies often suggest many questions for further 
study and we conclude this note by discussing some questions raised by this present study. 

First it must be noted that the powers of the poles at any particular location increase 
with the power of k. This is an indication that the location of the true singularity depends 
on k and it would be most desirable to find a perturbation theory directly for the location 
of the singularities. 

In terms of extending our results we remark that while it may be cumbersome to extend 
the results for the non-integrable case it is clear that for the integrable and superintegrable 
cases a great deal more can be computed for the perturbation expansion of the correlation 
functions. It is to be hoped that further terms in these expansions will produce conjectures 
such as ( |3.1| ) for the full correlation functions. 

Finally there is the question of whether or not the increased analyticity which results 
from the the cancelation of poles can be taken as a sufficient condition for integrability 
which can be considered as part of problem of the general relation between analyticity and 
integrability. 

The construction of models which satisfy the definition of integrability ( |1.5| ) in terms 
of commuting transfer matrices can be viewed as essentially a problem in algebra. Indeed 
one of the major developments in algebra in the last 15 years is quantum groups which 



were invented specifically in order to find classes of solutions to (|1.5| ). Analysis, on the 
other hand seems to have only a marginal relation to ( |1.5|) and, in the larger scheme of 
mathematics, algebra and analysis are almost universally treated as separate disciplines. 

A corresponding separation of disciplines has been traditional in physics where sta- 
tistical mechanics has relied heavily on algebra while field theory has relied on analysis. 
Thus, for example, while Onsager's solution of the Ising model is totally algebraic the 
renormalization theory of quantum electrodynamics is almost exclusively concerned with 
subtractions of infinities in integrals and seems to be totally concerned with analysis. 



11 



Thus 50 years ago integrability (algebra) and analyticity (analysis) were seen to be 
totally separate. However, in the 50 years that have followed these two totally separate 
computations there has been a remarkable merging of field theory and statistical mechanics 
which can be viewed as either the introduction of algebra into field theory or the introduc- 
tion of analysis into statistical mechanics. Thus, as an example, at the present day we have 
the isomorphism in two dimensions of solvable statistical mechanical models (solved by al- 
gebra) with conformal quantum field theory (where analysis appears in the very statement 
of the subject) . 

If we take this merging of analyticity and integrability seriously, however, then we 
are forced to conclude that the recent four dimensional studies on = 2 Yang Mills the- 
ory which have total reliance on analyticity properties should have an isomorphic 
counterpart whose solution relies on algebra. The discovery of such an algebraic counter- 
part (or the demonstration that it does not exist) to the work of will constitute an 
important advance in our understanding of the relation of integrability and analyticity. 

Acknowledgments. The authors wish to thank Prof. Gert Almkvist, Prof. G. von Gehlen, 
Prof. A. Guttmann and Dr. A. Honecker for insightful discussions. This work was partially 
supported by NSF grant DMR 9404747. 



Appendix A. Proof of Trigonometric Identity ( [4.12| ) 



To prove ( |4.12| ) we begin by using the symmetries of f{N,i,k) and the cotangent 



function to write 

I N-l I I JV-1 I 

S{N)=A J2 E Wj,^)cot^cot^+4 P(iV,j,j)cot2^. (Al) 

j=l k=j+l j=l 

We recall that P{N,j, k) was defined by (|4.14|) in terms of a function Q{N, x) satisfying 
the symmetry property ( [4.15| ) and then write Q{N,x) so that the symmetry property 
automatically holds 

Q{N,x) =r{x) -r{l- x) (A.2) 
12 



where r{x) is an arbitrary function whose dependence has been suppressed for brevity. 
Using this definition we expand ( |A.1|) as 



7V-J-1 r 



j=i k=j+i 



+4 E 

[^\ N-j-l 

j=i k=j+i >- 



k+j 
N 

+ r(0) + 2r 

N-k-j 
N 



nj nk 
cot -cot - 



N-j 
N 



cot 



2 ^ 



iV- A; 



cot -cot - 



^ E 



which may be rewritten as 



N-2j 
N 



+ r(l) + 2r 



cot^ 



A^ 



a=0 



a=0 



where 



4E[=r ^cot 



2 TTJ 

N 



for a = 



with 



CiN, a) - <j ^Ci{N, a) + 4(^2 (A^, a) for 1 < a < [^J 
-C(A^, N-a) for aU a 



|_ 2 J . , a-1 

nj 7i[a — j) o J- ^ 7ra 



Ci(A^,a)=2 ^ cot cot 



A^ 



\ ^ 71 J 71 U 

2 > cot — cot h ea cot 

^ A^ A^ 



Tra 
2iV 



and 



a-l 



C2(A^, a) = - 2 ^ 77 



— 1-2 > cot — cot 



7ra 



— epj-a cot^ 2 cot^ 

2A^ 



Here 



A^ 

Tra 
iV' 

1 if a is even 
if a is odd. 



A^ 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 
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The sum C{N, 0) may be rewritten as 



N-l 

C{N, 0) = 2 ^ cot 



2 ^ 

N 



which is recognized as the Dedekind sum and hence 

3 

To evaluate Ci{N, a) we first write 

E 



^ / AT \ 7r(a — j) ^ n(a — j) nj 

Ci{N, a)= cot cot ^ ^ + 2^ cot ^ „ ^ cot 



N 



N 



2 7ra ■^-^ /I J 
+ ea cot ^-22^ cot — cot 



a-l 



2N 



nj 7ra 
— cot — 
N N 



a — l . I -x a — 1 

ETTj 7r(a — TTj 7ra 
cot — cot 2 > cot — cot 



a-l 



$:cot^cot"(^'+^-"^ 



iV 



iV 



- 2 cot -^^T — - cot 



N 



Then let u — exp (^) and write the first summand as 

Tvj Tv(j + N -a) uj^ + 1 uj^+^-'' + 1 

— cot cot = : ■ —TT 

uj^ -1 - 1 



A^ 



1 + 
1 + 
1 + 
1 + 
+ 



UJ3 -1 



1 + 



+ 



2 



+ 4 — 



+ 



^3 - 1 - 1 

2 

- - 1) 

2 



+ 4 



1 



N- 



1 



^w-?' - 1 - 1 / w^-" 

(a;^-"-l + 2) 



AT-a 



- 1 - 2w 



Af-a^ 



1 + 2- 



JV-a 



+ 1 / 1 



1 -1 - 1 
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The sum partially telescopes 

o-l 

E n J 
cot — cot 
N 



TTJ Tl{j + N — a) 

w^-" + 1 



(a - 1) + 2 



- 1 

1 



1 1 

+ ^ r + --- + 



1 



1 



N-a 



- 1 



(a-1) + 2 



u^-" + 1 



- 1 



,Ar-a 



(a - 1) + 2cot 



7r(A^-a) 



uj^ + 1 uj'^ + 1 u^-" + 1 

+ r H + 



UJ^ - I ' - 1 

N-a 



N-a 



- 1 



E 



cot 



It follows that 



Ci{N, a) = a - 1. 
A parallel computation shows that 



N-a-l 



C2(iV,a)= 



" — 2 cot — > cot — 



cot 



iV 



and therefore for a ^ 



Ci{N,N-a) = -(N-a-l) 
C{N,a) = -4{N-2a). 



(A.13) 



(A.14) 



(A.15) 



(A.16) 



Appendix B. Additional Identities Related to ( |4.12| ) 

The cotangents on the left hand side of (|4.11[) may be expanded in powers of iv = 
exp (^) but the roots of unity all cancel in the sum. We expect that similar cancelations 
will occur when the sum rule for the superintegrable correlation functions is applied to 
orders higher than four. Thus we have investigated the generalization of (|4.11| ) 



where 



N-lN-l . , 

S{N) = hp{N,i, k) cot — cot — = polynomial in 

j=i k=i 



( P{N, j, k) if 1 < i < A; < A^/2 

hp{N , i,k) = < hp{N,k,j) in general 

y—hp{N,N — j,k) in general 



(B.l) 



(B.2) 
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and where P{N, j, k) is a polynomial in A^, j and k of total degree d. In the case where is 
prime, the only way for the roots of unity to vanish is for the coefficients of cu, uj^, . . . , uj^~'^ 
all to be equal. This implies a system of equations for P{N,j, k) which is overdetermined 
when is sufficiently large. Remarkably, there are many solutions to this system and 
these solutions hold for all values of tested. We conjecture that they are true for all A^. 

Up to the largest degree tested, d = 14, the number of solutions to the overdetermined 
system is equal to the number of solutions to 

d = 3f + 2n (B.3) 

with / a positive integer and n a non-negative integer. Thus we make the further conjecture 
that there are infinitely many solutions which can be organized into families, / = 1,2,3,..., 
whose members, labeled by the index n, have degree given by ( [B.3| ). Although there is 
immense freedom in combining solutions, it has proven to be possible to write the solutions 
in such a form that general formulas, given below, can be guessed in the cases / = 1, 2. 
The solutions, Pj^j^(A^, j, k), that we have found for / = 3,4 are shown in Table 8 and the 
corresponding sums, Sf^n{N,j, /c), are shown in Table 9. 

The first family, / = 1, comprises polynomials of the form (|4.14|) . If in equa- 
tion ( [4.12|) we make the choice, Q{N,x) = N'^qd^x) with qd{x) a polynomial of de- 



gree (i, then S{N) is a polynomial in A^ of degree d -\- 2. The symmetry property on 
Q{N,x), ( |4.15| ), requires d to be odd and d = 1 implies P{N,j,k) = 0. Thus these 
solutions have the degree stated in ( p.3|) . 



The second family of solutions is given by 



j \ 2n+3-2h , ^ 

\hB2h+3 1 ]y ) ^" 



h=0 

■ \ 2n+2-2h / ^ 

where, letting g{n) = 4"+^ — 1, 

/2n + 3^^ 1 g{h)g{n-h) 



(B.4) 



2h J {2h + l){2h + 2){2h + 3) g{n) 
2n + 3\ 1 g{h)4''-^ 



(B.5) 



2h + lJ {2h + 2){2h + 3) g{n) 
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and where Bn{x) and En{x) are the BernouUi and Euler polynomials, 
in ( p.4| ) may be performed using the identities 



The sum 



Br,{x + h) = J2Cf)Bk{x)h 

k=0 ^ ^ 
I — n \ / 



(B.6) 



to obtain 



2n 



{2n + 6)(2n + 5)(2n + 4)(4"+i - 1) 



2n+6 



k + j 
N 



(6 ■ + l)S2n+6 ( ^ ) + (4"+' - l)B 



2 . 4^+2 / 5 



2n+6 



2 . 42^^+4 ( S 



2n+6 



2-4-+2(S2n+6 



k+j/2 

N 
k + j 
2N 

k + 2j 
2N 



2n+6 1 jY ) + -S2n+6 . jy 



k-2j 



B2n+6 ( M— I I + 



+ B2n+6 
+ B2n+6 



N 
k-j 



2N 

k-2j 
2N 



V\B 



2n+6 



2N 



(B.7) 

This expression differs from ( |4.14| ) in that the Bernoulli polynomials cannot be replaced 
by arbitrary polynomials satisfying g(l — x) = —q{x). A general form has not been found 
for the corresponding sums, S2,n{N). The first five are shown in Table 7. 
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Table 1 The first nine orders of perturbation in k for N = 3 for the ground state 
energy in the generic case. We use the variable C3 = cos (0/3). 



(3) _ (3) _ -2 

(3) _ -COS0 

e(3) = 1 + ^ ! 

^ VI62C33 8IC3 81(-1 + 2C3) 81(1 + 2C3)/ 



g(3) _ 



/- / — COS0 4cos0 8cos(/) 8cos0 \\ 

V108C3^ ~ 81C^ ^ 81 (-I + 2C3) ~ 81 (1 + 2 C3) 7 y 



J3) - ,/3 ( -2 1 2(l-54cos^<^) 

\729(-l + 2C3)' 729(-l + 2C3)' 2187 (-1 + 2 C3) 



1 2 (1 - 54 

+ — — — -o + ^ 



cos^ 



729(1 + 2C3)^ 729(1 + 2C3)^ 2187 (1 + 2 C3) 
ll + 24cos^^ 5 + 54cos^0 2 (-1 + 54cos^ ^) 
^ 11664 C3^ ^ 4374 C3^ ^ 2187 C3 

(3) _ I COS 263 cos 4> 1022 cos (t> 



^972C3^ 26244 C3^ 19683 C3^ 
COS0 64 cos 119 cos 49 cos ^ 



e 



2187(-l + 2C3r 6561(-l + 2C3f 39366 (-1 + 2 03)^ 486 (-l + 2C3) 
cos^ 64 COS 119 cos ^ 49 cos ^ \ 

^ 2187(1 + 2(73)'' ~ 6561 (1 + 2 C3)^ ~ 39366(1 + 2(73)^ ^ 486 (1 + 2(^3) J 

(3) _ ^ / -1 -1-2 cos^ 4> _|_ 55 - 888 cos^ ^ 



2187 (-1 + 2(73)^ 2187 (-1 + 2(73)^ 78732 (-1 + 2(73)^ 



-451 + 1384cos^0 2759 + 15264 cos^(?!) 1 1 + 2 cos^ (p 

^ 157464 (-1 + 2 Caf ^ 78732 (-1 + 2 (73) 2187 (1 + 2 Csf ^ 2187 (1 + 2 Cs)"" 

^3 



55- 888cos2(* 451 - 1384 cos^c^ 2759 + 15264 008^(7!) 1715 C3 



+ — — — — ^ + + 



78732 (1 + 2 Caf 157464 (1 + 2 Caf 78732 (1 + 2 (73) 2187 (-3 + 28 Ca^) 

-381 - 1936 cos^ ^ -801 - 6208 cos^ ^ -709 - 7880 cos^ 4> -277 - 7632 cos^ <? 
+ . + + + 



1679616(73^ 629856(73^ 157464(73^ 39366(73 
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(3) nz I 5 COS d) 17 COS 4> 61 COS 4) + 432 cos'^ 



26244 (-1 + 203)" 8748 (-1 + 2 Ca)-" 944784 (-1 + 2 C; 



3J 

-11257cos(^ + 6048cos^(^ 105767 cos - 24912 cos^ ^ -34607 cos ^ - 59184 cos^ 
^ 1417176 (-1 + 2 Caf ^ 5668704 (-1 + 2 Cs)^ ^ 629856 (-I + 2C3) 

5cos^ _ 17 cos ^ 61cos^ + 432cos^(^ 11257 cos ^ - 6048 cos^ 

26244(1 + 2(73)^ ~ 8748(1 + 2(73)^ 944784(1 + 2(^3)'' 1417176(1 + 2(^3) 

105767 cos 0- 24912 cos^^ 34607 cos + 59184 cos^ ^ 24010 cose 



3 



5668704 (1 + 2 (^3)^ 629856 (1 + 2 C3) 19683 (-3 + 28 Cs") 

133cos^+ 192cos^(^ 13 (ll cos + 52 cos^ 0) -575 cos ^ + 3696 cos^ 



+ 



839808 314928 C3'' 314928 (^3' 

-9007cos(^ + 17424 cos^. 



354294 (73 ' 
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Table 2 The first eight orders of perturbation in A; for AT = 4 for the ground state 
energy in the generic case. We use the variable C4 = -\/2cos(0/2). 



(4) _ cos</) COS0 cost;!) 



4C4 4(1 + C4)' 4(1 + C4) 



^(4) _ -1 1 5-2 cos^ 1 8 + 9 cos^ cf) 

~512C7^16C7^ 160 ^ 12 (1 + Cif ^ 72 (1 + C4)' 
_^2O + 27cos^0 875 



216(1 + C4) 432 (-l + 5C4) 

(4) _ 3cos^ 19 cos ^ 149 cos 61 cos ^ 55 cos 

" 128C4^ ^ 64C4' ^ 128 C4 ^ 576 (1 + C4)'' ^ 864 (1 + C4)^ 

cos 4) 277 cos 1225 cos ^ 89825 cos, 



,(4) 



1152 (I + C4)' 31104 (I + C4) 648(-l + 5C4)^ 15552 (-1 + 5 C4) 

2025 1 3 -7 + 52cos2<^ 3 (-23 + 48cos^ (^) 



+ . + 



8192(-1-3C4) 8192C4^ 512C4^ 1024C4^ 512C4 



2 



-303 + 812cos^(;i!) 155-44 cos^ cj) -33- 394cos^^ -253 - 3911 cos^ </> 



1024 C4 3456 (1 + C4)^ 3456(1 + ^4)^ 20736 (I + C4) 

-15832 - 72167 cos^t^ 75203 - 236166 cos^ 8575 (5 + 4 cos^ (^) 



^ 373248 (1 + C^f ^ 1492992 (1 + C4) 15552 (-1 + 5 C4)^ 

175 (8345 - 12704 cos^ $) 25 (727595 - 758184 cos^ 9(1 + Ca) 



746496(-l + 5C4)^ 5971968 (-1 + 5 C4) 16 (l - 2 C4 - 2 



(4) _ 5cos(/) 31 cost;/) 1301 cos — 96 cos 309 cos </) + 20 cos^ 
^ ~ 4096C?" ~ T024C^ 4096C? 256 

3855 cos ^ + 2944 cos^ ^ 101 cos cj) + 960 cos^ ^ 13385 cos ^ + 5904 cos^ ^ 



4096 C4 41472 (1 + C^f 124416 (1 + C4)^ 

124655 cos ^ + 83952 cos^ 4> _ 96257 cos ^ + 14169 cos^ _ -4392121 cos ^ - 488304 cos^ ^ 

1492992 (1 + CiY 559872 (1 + df 26873856 (1 + df 

11116697cos^- 1116063 cos^^ 1539cos(^ _ 25839cos(^ 

20155392 (1 + C4) ^ 16384 (1 + 3 C4)' ' 32768 (I + 3C4) 

1860775 cos ^ 1225 (4103 cos ^ + 2800 cos^ ^) 



i3 



186624 (-1 + 5 C4) 746496 (-1 + 5 C4 

25 (-47222099 cos + 621600 cos^ 25 (432042539 cos ^ + 85572000 cos^ ^) 

35831808 (-1 + 5 C4)^ ^ 644972544 (-I + 5C4) 

51 (2 cos + cos C4) 
16 (-1 + 2C4 + 2C4^) 
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81 (199 + 32 008^0) 81 (-803 + 7760 cos^^) 81 (-20863 + 102076 cos^ 
~ 131072 (-1 - 3 ^ 1048576 (-1 - 3 C4)^ ^ 13631488 (-I-3C4) 
25 9 -151 - 158cos^^ -519 - 3008 cos^ 



+ 



2097152 C4^ 16384 C4^ 16384 C4^ 16384 C4 

359 - 15946 cos^ <^ - 32 CDs'* 4> 



16384 C4' 



7923 - 43936 cos-^ (p - 2048 cos^ (p 76123 + 14074 cos-^ cp - 10816 cos^ (j) 
^ 16384 C4 2 ^ 16384 C4 

-53297 + 31488 cos^ 4> - 4608 cos"^ 4> -555061 - 193632 cos^ ^ + 138240 cos^ . 



+ 
+ 



1990656 (1 + C4)^ 23887872 (1 + 

-1933609 + 12161664 cos^ ^ + 2257920 cos^ ^ 

95551488 (1 + 
85101083 - 62007936 cos^ ^ + 179776512 CDs'* 



3439853568 (1 + C4 



,4 



-2291995019 + 25564353792 cos^ ^ + 3332192256 cos^ 4> 
41278242816 (1 + CUf 

9085861 571 - 4' 

+ 



9085861571 - 47135217152 cos^ 6 + 5619560448 cos"^ . 



55037657088 (1 + C4? 



_^ -917428110137 + 7122467770368 cos^ (p + 222666719232 cos"^ cp 420175 (2 + 3 cos 



1981355655168 (1 + C4) 69984 (- 1 + 5 C4 

8575 (-3751 + 21936 cos^ <^) 175 (161445421 + 287088144 cos^ ^ - 35280000 cos"^ 
"l~ , \ _ _ > /I "l~ 



,5 



+ 



+ 



26873856 (-1 + 5 C4)'* 2579890176 (-1 + 5 Cif 

175 (-835186517 + 7197054312 cos^ 4> + 37920000 cos^ $) 

10319560704 (-1 + 5 C4)^ 
125 (1145946953233 + 5452562686272 cos^ + 282187776000 cos'' 0) 
12878811758592 (-I + 5C4) 
29176875 9 (4 + 3 C4) 3 (81 + 46 C4) 



+ 



33554432 (1 + 5 C4) 128 (l - 2 C4 - 2 C^)' 512 (l - 2 C4 - 2 C^)' 

- 19777 + 197496 cos^ ^ - 22351 C4 + 197496 cos^ ^ C4 27 (48104891 + 84022143 C4) 



6656 (1 - 2 C4 - 2 C4^) 16777216 (-1 + 8 C4 + 17 C4^) 
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Table 3 The first 5 orders of perturbation in k for N = 5 for the ground state energy 
in the generic case. For the fourth and fifth orders only the terms proportional to cos^ 
and cos^ respectively are given. We use the variable C5 = '^°^^/5'' • 



e(^) = -(5- V5)C| + (7-y5)C5, ef^=0 

e(^) = (5 + V5) ( -(^ + ^) I 9 + 75 3 + 75 {-3 + VE) 

{^^^^) \5Q(^_l + C,)^ 5OC5 50 (1 + Cs) ^ 50 (-2 + Cs^) 



^5 + / - ( (5 + 3 VS) COS 0) (35 + 13 V5) cos(^ V2 (25 + 6 V5) cos ^ 
~ -5 + 3^5 V 125 V2 {-I + C5f ^ 125^2 (-1 + C5) 125C? 
(5 + 3^5) cos ^ (35 + 13^5) cos ^ 2^ (-5 + 2 VS) cos ^2 (5 + 4 ^5) cos , 



125^2(1 + ^5)' 125^2(1 + ^5) 125(-2 + C52)' 125 (-2 + ^5") 

(5) _ 2 7 ( 110 + 50^/5 _ 4 (15 + 7^5) 260 + 120v/5 630 + 290 VS 
~ V3125(-l + C5)' ~ 625 (-I + C5) + 3125 Cs^ ^ 3125 C5 

110 + 50 4(l5 + 7^/5) (-40 - 20 Vs) C5 (30 + 10 Vs) C5 



3125 (I + C5)' 625 (1 + ^5) 3125(-2 + C52)' 3125 (-2 + 

(5) _ ( -52 V1O-2V5 _ 116 VlO-2V5 \ 3 - 

3125 3125^ j V4(-l + C5)' 4(-l + C5) 

13 1 1 1 3 

+ + + — — o + . . ^ s + 



4C5^ 4C5' 4(1 + C5)' 4(1 + C5) 4(_2 + C5')' 4(-2 + C52) 
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Table 4 The first 6 orders of perturbation in k for N = 6 for the ground state energy 
in the generic case. For the sixth order only the term proportional to cos^ is given. We 
use the variable Cq = -\/3cos(0/3). 



J6) ^ -2 _ A + 64 3 

' 15(-1 + C6) Ce 15(3 + 2C6) 4(-9 + 8C6^) 

(6) —cose; 91cos(p cos4> 8cos0 

~ 75 V3(-l + C6)^ 375V3 (-l + Ce) VSCe^ SVSCe 

484 cos 132^3 cos 32cos(^C6 

~ 75V3(3 + 2C6)^ ~ 125 (3 + 2C6) ~ 3^3 {-9 + 8Ce'') 

(6) 1 667 + 360cos^(^ 216617 + 527040 cos^ 

~ 2250 (-1 + Cef ^ 15000 (-1 + Cef ^ 450000 (-1 + Ce) 

11 2 (-39 + 50 cos^ (^) ^ -3691 + 5000 cos^ 
^ 30C7 ^ 2257V 6750 Ce 

2 (-2381 + 3060 cos^ ^) -1620919 + 579960 cos^ ^ _ 708766379 + 423636480 cos^ ^ 
^ 3375 (3 + 2C6)^ ^ 911250 (3 + 2 Ce)^ 492075000 (3 + 2C6) 

1294139 2 (3 + Ce) 4 (9 + 8C6) 

~ 243000 (-15 + 14C6) " 99 (-4 + C6 + 2C6') ~ 297 (-3 + 4C6') 

27 4 (-9 + 16C6) _ 44 (-11 + 36 cos^ + 8 Ce) 

~ 64(-9 + 8C62)' 9(-9 + 8C62)' sU^d + SC^) 

^ 9878103 - 19929087 Ce + 7810406 Ce^ 
314928 (27 - 36 Ce - 33 Ce^ + 38 Ce^) 



23 



(6) _ COS (j) 193 cos (p cos (f) (-242459 + 194400 cos^ (f) 

" 1000 V3(-l + Ce)^ 10000 ^ (-1 + C6)^ 600000 ^ (-1 + Ce) 

5427^3 cos ^ 126279 cos 3 ^3 cos (-185083 + 11250 cos^ (^) 
^ 160 Ce^ 800 Ce^ 2000 Ce 

81 73 cos ^ (-65789 + 5400 cos^^) 9y3cos<^ (-842051 + 361575 cos^ 
1000 (3 + 2 Cef 2500 (3 + 2 Cef 

cos^ (-2735083651 + 191181600 cos^ (f) 8349 a/3 cos 

^ 800000^ (3 + 2C6) ^ 32000 (-15 + 14^6) 

_ 469 cos ^ (-302 + 255 Ce) _^ VScosc^ (-229 + 264 Ce) 

186 73 (-4 + Ce + 2Ce') 8 (-3 + 4Ce') 

_^ 3645 73 cos (-123 + 116 Ce) _ 27 ^3 cos (-83589 + 79246 Ce) 
1024 (-9 + 8 Ce^)^ 512 (-9 + 8 Ce^)^ 

9 73 cos ^ (8025965 - 2928384 cos^ - 8380680 Ce + 2761728 cos^ 4> Ce) 

31744 (-9 + 8Ce^) 
cos ^ (22455 - 53739 Ce + 29648 Ce^) 



4 73 (27 - 36 Ce - 33 Ce^ + 38 Ce^) 



27 1674 1 



^ = cos d) \ + / . . ^ s + 



625(-l + C6)^ 3125 (-1 + Ce) 27C6^ 81 Ce 
16 64 13312 



75(3 + 2Ce) 125(3 + 2Ce) 16875 (3 + 2 Ce)' 
258688 256 



253125 (3 + 2 Ce) 27 (-9 + 8 Ce^) 



+ 
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Table 5 

The results of specializing the generic results of the perturbation theory for the ground 
state energy for = 3, 4, 5 and 6 to the integrable manifold (|1.4|) where cos (j) = k cos cj). 



^(3) _ (3 + 4C3^) _ ^,{3 + 16Cs^ 

^ \/3 9V3C3 324V3C33 

,6 27 + 108C3^ -72(:73^ + 800C3^ , „,,g 
-k ^ — h 0[k 

52488 V3C3'^ 



4 8 2 (1 + C4] 



' 64 + 512 + 32(1 + C,f ^ 64(1 + ^ ^ ^^''^ 



A5) 



(5 - V5)Cl + (7 - V5)C5 



+k' I ^ + —r^ + 



1 2 (-3 + ^) 



,5 (-5 + 2^) (-I + C5) 5(-5 + 2^)C5 5 (-5 + 2 ys) (1 + Cs) 
C5 (18-8V5-7C5^ + 3V5C5^) 



5 (-5 + 2^5) 



+ 0(fc^) 



... .3 - ICe - Id - e^yi + ±C. + 5/36 + ± - + O,.' 
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Table 6 The first six orders of perturbation theory in k for N = 3 for the correlation 
function (ZqZ^) = X]n=o ^"^n generic case. We use the notation C3 = cos 0/3 and 

S3 = sin 0/3. For cf only the coefficient of 5r,i and for Cg only the coefiicient of cos^ 0(^r,i 
is given. 



Co = 1, Cf = 0, C2 = ^((5,.,o - 1), ^3 = Y^('5r,0 - 1) 

2 2 2 2 ^ \ /i £ \ 

^"^^ '^432C|^ 81^~ 27(1-2(73)2 ~ 81(1 - 2C3) ~ 27(2C3 + l)^ ~ 81(2C3 + 1)^*^ ~ 
/ 5 1 2iS3 

M32C4 ^ 81C| ^ 27C3 

1 - 2iS3 7 - 12i53 1 + 2ig3 7 + 12ig3 w 

54(2C3 + 1)2 162(2C3 + 1) 54(2C3 - 1)^ 162(2C3 - 1)^ 

z_ T/^ 1,2, 2ig3 

Cg - cos ^Q^scl 81CI ^ 8IC3 + 27C| 

1 + 2iS3 2 + 18^53 -1 + 2x53 -2 + 18x53 x 

27(2C3 - 1)2 ~ 81(2C3 - 1) ^ 27(2C3 + 1)2 ^ 2C3 + 1 ^'^ + ' ' ' 



2Y(2G3 - J 

2 7/ -1 1 1 i53 , 4i53 



t t ^ + ^^ + 2!^ 

1296Cf 108C| 54C2 54C| 27C3 

1 - 2iS3 1 + SiSs 1 + 2iS3 

54(2^3 + 1)2 54(2(^3 + 1) ^ 54(2C3 - 1)2 ^ 54(2C3 - 1) ^ ^ 
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Table 7 The first five sums, ^2,n(A^), from ( [gT| ) with P2,n(A^, J, k) given by (|B]|). 

52,o(A^) = -N^{N^ - 1){N^ - 4)(53iV^ + 13)/(16 • 27 • 5 • 7) 
52,i(7V) = N^{N^ - 1){N^ - 4)(821Ar^ + 2807V^ + 69)/(8 • 81 • 25 • 7) 

S2,2{N) = -N^{N^ - l)(iV^ - 4)(308747V^ + lllSOiV^ + 3501iV^ + 845)/(32 • 81 • 5 • 7 • 11) 
S2,3iN) = N^{N^ - 1){N^ - 4)(103801414Ar^ + 37967420A^^ + 12484059iV'^ + 3718690^"^+ 

882407)7(8 • 27 • 125 • 7 • 11 • 13 • 17) 
S2,4{N) = -N^{N^ - 1){N^ - 4)(10128697897V^° + 371602425iV^ + 123446295iV'^ + 

38256875/^^"* + 11003386A^^ + 2579010)/(32 ■ 81 • 5 ■ 7 ■ 11 • 13 • 31) 
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Table 8 The first three members of the third family and the first two members of 
the fourth family of solutions to (p.l| ) expressed in terms of the Bernoulli polynomials and 
using the notation J = j /N, K = k/N. 



P3,o(iV, j, k) = N^{J^{K - 2K^ + K^) - 2J^{K - K^) + J^{K + l/?,K^)) 
P3,i{N,j,k) = N^^{l/W^{K - 6K^ +9K^ - AK"^) + 3/2J^{K^ - K^)- 

l/AJ^{QK - 9K^ + 7K^) + 2f{K - K^) - 1/4J^(3K + K^)) 
P3,2iN,j,k) = N^\J\5/9B9{K) + Bj{K)) + j\l/l2B-j{K) - l/mBs{K)) + 

j\2B-j{K) + Al/24B5iK) - l/SBsiK)) + f?,/AB^{K) + 

J\A/2,B^{K) + Al/72Bs{K) - l/16fii(if)) + 7^5/9^3^+ 

J^°(l/953(A') + l/9Si(i^))) 
P4,o{N,j,k) = N^^{J^{K -3K^ + 3K^ - K'^) - 3J^{K - 2K^ + K^) + 

f{3K - 2K^ - K^) - J^{K + K^)) 
P4,i{N,j,k) = N^^{j\b/9B9{K) + Br{K)) + j\2/3B7{K) + l/2B^{K) - l/lSBs{K)) + 

j\U/9Br{K) + 17/9B5{K)) + J^{4/3B5{K) + 1/2B3{K) - l/36Bi{K)) + 

J\b/9B^{K) + 20/27^3 + J^\2/9B3{K) + l/9Bi{K))) 
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Table 9 The sums ( |B.1|) corresponding to the polynomials in Table 8. 

S3,q{N) = N{N^ - 1){N^ - 4)(iV^ + 5N^ + BiN^ - 80) • 32/(27 •25-7 
53,1 (AT) = N{N'^ - 1){N^ - A){9N^ + 4.5N^ + A6N'^ + 12920iV^- 

31680) • 8/(27-25 • 7 - 11 • 13) 
S3,2{N) = -N{N^ - 1){N^ - 4)(198A^^° + 990A^^ + 7980A^^+ 

581057V'' + 15362A^^ - 786240)/(64 • 243 • 25 • 49 • 11 • 13) 
S4,o{N) = N'^{N^ - l)(iV^ - 4)(27V* + lOTV*^ + 42A^'' + 2315A^^- 

9614) -256/(27 -25 -49 -11 -13) 
5'4,i(A'') = -A^^(A^^ - l)(iV^ - 4)(97V^° + 45A^® + 1897V'^+ 

79157V'' + 11792A^^ - 240240)/(32 • 243 • 25 • 49 • 11 • 13) 
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